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We present a method for evaluating divergent non-Borel-summable series by an analytic continua- 
tion of variational perturbation theory. We demonstrate the power of the method by an application 
to the exactly known partition function of the anharmonic oscillator in zero space-time dimensions. 
In one space-time dimension we derive the imaginary part of the ground state energy of the an- 
harmonic oscillator for all negative values of the coupling constant g, including the non-analytic 
tunneling regime at small —g. As a highlight of the theory we retrieve the divergent perturbation 
expansion from the action of the critical bubble and the contribution of the higher loop fluctuations 
around the bubble. 



I. INTRODUCTION 

None of the presently known resummation schemes 0, Q is able to deal with non-Borel-summable series. Such 
series arise in the theoretical description of many important physical phenomena, in particular tunneling processes. 
In the path integral, these are dominated by non-perturbative contributions coming from nontrivial classical solutions 
called critical bubbles 0,0 or bounces and fluctuations around these. 

Any Borel-summable series becomes non-Borel-summable if the expansion parameter, usually some coupling con- 
stant g, is continued to negative values. Here we show that non-Borel-summable series can be evaluated with any 
desired accuracy by an analytic continuation of variational perturbation theory 0j in the complex g-plane. This 
implies that variational perturbation theory can give us information on non-perturbative properties of the theory. 

Variational perturbation theory has a long history @, 0, IE El It is based on the introduction of a dummy 
variational parameter Q on which the full perturbation expansion does not depend, while the truncated expansion 
does. An optimal f2 is selected by the principle of minimal sensitivity |lOj|. requiring the quantity of interest to be 
stationary as function of the variational parameter. The optimal Q is usually taken from a zero of the derivative with 
respect to tt. If the first derivative has no zero, a zero of the second derivative is chosen. For Borel-summable series, 
these zeros are always real, in contrast to statements in the literature 0, 0, fiH. Il4| which have proposed the use of 
complex zeros. Complex zeros produce in general wrong results for Borel-summable series, as was recently shown in 
Ref. 

The purpose of this paper is to show that there does exist a wide range of applications of complex zeros if one wants to 
resum non-Borel-summable series, which have so far remained intractable. These arise typically in tunneling problems, 
and we shall see that variational perturbation theory provides us with an efficient method for evaluating these series, 
rendering their real and imaginary parts with any desirable accuracy, if only enough perturbation coefficients are 
available. An important problem which had to be solved is the specification of the proper choice the optimal zero 
from the many possible candidates existing in higher orders. A non-Borel-summable series is associated with a function 
which has an essential singularity at the origin in the complex (/-plane, which is the starting point of a left-hand cut. 
Near the tip of the cut, the imaginary part of the function approaches zero rapidly like exp(— for g — > 0—. If 
the variational approximation is plotted against g with an enlargement factor exp(cv/|(7|), oscillations become visible 
near g = 0. The choice of the optimal complex zeros of the derivative with respect to the variational parameter is 
fixed by the requirement of obtaining, in each order, the least oscillating imaginary part when approaching the tip of 
the cut. We may call this selection rule the principle of minimal sensitivity and oscillations. 

In Section ^ we shall explain and test the new principle on the exactly known partition function Z{g) of the 
anharmonic oscillator in zero space time dimensions. In Section IIIII we apply the method to the critical-bubble 
regime of small — g of the anharmonic oscillator and find the action of the critical bubble and the corrections caused 
by the fluctuations around it. In Section I1VI we present yet another method of calculating the properties of the 
critical-bubble regime. This method is restricted to quantum mechanical systems. Its results for the anharmonic 
oscillator give more evidence for the correctness of the general method of Sections [D] and IIIII 
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II. TEST OF RESUMMATION OF NON-BOREL-SUMMABLE EXPANSIONS 



The partition function Z(g) of the anharmonic oscillator in zero space-time dimensions is 



z{g) 



^ v i-x 2 /2-g x 4 /4) dx = CXP ^ 9 ) -K 1/4 (l/8g) 



(2.1) {FOKKER} 



where K v (z) is the modified Bessel function. For small g, the function Z(g) has a divergent Taylor series expansion, 
to be called weak- coupling expansion: 



32*00= £<w. with ai = (-i)' r(2 !, + / 2) . 

For g < 0, this is non-Borel-summable. For large |<?| there exists a convergent strong- coupling expansion: 

z£Ug) = 9 ~ l/i £ 6, <r i/2 , with bl = (-iy r(V 2 2 , + ' /4) . 



(2.2) {FP-WEAK] 



(2.3) {FP-STR0I> 



(=0 



As is obvious from the integral representation Ij2.1|l . Z(g) obeys the second-order differential equation 

l6g 2 Z"(g) + 4(1 + 8g)Z'(g) + 3Z(g) = 0, (2.4) 

which has two independent solutions. One of them is Z(g), which is finite for g > with Z(0) = oq. The weak- 
coupling coefficients a; in (|2.2() can be obtained by inserting into l|2.4[) the Taylor series and comparing coefficients. 
The result is the recursion relation 



{DGL} 



ai+i 



161(1 + 1) 



4(2 + 1) 



■ai. 



(2.5) {REC} 



A similar recursion relation can be derived for the strong-coupling coefficients 6/ in Eq. I|2.3[) . We observe that the 
two independent solutions Z(g) of 1|2.4(1 behave like Z(g) cx g a for g — > oo with the powers a = —1/4 and —3/4. The 
function (|2.1|) has a = —1/4. It is convenient to remove the leading power from Z(g) and define a function £(x) such 
that Z(g) = g- 1 ^ C{g~ 1/2 ). The Taylor coefficients of ((x) are the strong-coupling coefficients bi in Eq. (|2.3I) . The 
function C( x ) satisfies the differential equation and initial conditions: 



4C"(ar) - 2a<'(a;) - ((x) = 0, with C(0) = b and C'(0) = h. 
The Taylor coefficients 6; of ((x) satisfy the recursion relation 

21 + 1 



5i+2 



-b,. 



4(l + l)(Z + 2) 

Analytic continuation of Z (<?) around g = oo to the left-hand cut gives: 

Z{-9) = (-.9)- 1/4 C((-.9)- 1/2 ) 

= (-.9)- 1/4 E^(-5)-' /2 exp 



(2.6) {DGL2} 



(2.7) {REC2} 



-^(21 + 1) 



for g > 0, 



so that we find an imaginary part 

oo 

ImZ(-g) = -(4. 9 )- 1 /4^ &; (_ g )-V2 sin 

oo 

= -(4.g)- 1/4 

where 

/3 = &Q, /3l = &1, A+2 = - 



(7T 



(2Z + 1) 



-1/2 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



{CUT} 



{CUTi} 



2/ + 1 



4(1 + l)(l + 2) 



(2.12) {CUT2} 
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It is easy to show that 

oo 

J2&x l =CWexp(-x 2 /4), (2.13) {«} 



so that 



ImZ(-g) = — ^.<T 1/4 exp(-l/4g) ^ h g- l/2 . (2.14) {CUT-STRC 



i=0 

From this we may re-obtain the weak-coupling coefficients a; by means of the dispersion relation 



Z(g)=-~ r imZ{ ~ z) dz (2.15) {DISP} 

7T Jo z + g 

1 f 00 exp(-l/4z) z-j/ 2 -V 4 



dz. (2.16) 



Indeed, replacing l/(z + g) by J °° exp (— x(z + g)) dx, and expanding exp {—xg) into a power series, all integrals can 
be evaluated to yield: 



oo oo 

Z(g) =- E (-3)'r(Z + J'/ 2 + 1/4) ■ (2-17) {DISP-2} 



7T 

3=0 1=0 



Thus we find for the weak-coupling coefficients ai an expansion in terms of the strong-coupling coefficients 

— 



a/=- V 2 3 &j T(l+ j/2+ 1/4). (2.18) {DISP-3} 

IT ^ — ' 

Inserting bj from Eq. i|2.3[l . this becomes 

ai = yijl £ r 0V2 + l/4)r(i + j/2 + 1/4) = (-1)' r(2 !, + ^ /2) , (2.19) {DISP-3} 

coinciding with 12.21) . 

Variational perturbation theory is a well-established method for obtaining convergent strong-coupling expansions 
from divergent weak-coupling expansions in quantum-mechanical systems such as the anharmonic oscillator 0, ITr| as 
well as in quantum field theory 0] • We have seen in Eq. (|2.8f) , that the strong-coupling expansion can easily be 
continued analytically to negative g. This continuation can, however, be used for an evaluation only for sufficiently 
large \g\ where the strong-coupling expansion converges. In the tunneling regime near the tip of the left-hand cut, the 
expansion diverges. In this paper we shall find that an evaluation of the weak-coupling expansion according to the 
rules of variational perturbation theory continued into the complex plane gives extremely good results on the entire 
left-hand cut with a fast convergence even near the tip at g = 0. 

The -Lth variational approximation to Z{g) is given by (see 

L 

Z(£(g, Q) = tT £ Ci(tr), (2.20) {FP-VAR} 



3=0 



with 



where q — 2/lu = A, p = — 1 and 



a = Cl q - 2 (Q^ - l)/g, (2.21) {FP-s} 



1=0 
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To apply the principle of minimal sensitivity, the zeros of the derivative of Zvai • (<?, with respect to are needed. 
They are given by the zeros of the polynomials in a: 

P {L) (<7)=Y,ai(p-lq + 2l-2L)f^ P L lq ^ 2 \-<7) L - l =0, (2.23) {FP-DERI\ 

1=0 ^ ' 

since it can be shown [l8| that the derivative depends only on a: 

Consider in more detail the lowest non-trivial order with L = 1. From Eq. I|2.23|l we obtain 




FIG. 1: Result of the 1st- and 2nd-order calculation for the non-Borel-summable region of g < 0, where the function has a cut 
with non- vanishing imaginary part: imaginary (left) and real parts (right) of zi^l(g) (dashed curve) and z!^ T (g) (solid curve) 
are plotted against g and compared with the exact values of the partition function (dotted curve). The root of giving 
the optimal variational parameter Q has been chosen to reproduce the weak-coupling result near g — 0. 



a =-, corresponding to f2 = - [\ ± ^/T+ lOgj . (2.25) {Zero} 

In order to ensure that our method reproduces the weak-coupling result for small g, we have to take the positive 
sign in front of the square root. In Fig. ^ we have plotted Z^l(g) (dashed curve) and Z^ [ (g) (solid curve) and 
compared these with the exact result (doted curve) in the tunneling regime. The agreement is quite good even at 

these low orders Next we study the behavior of zi^Hg) to higher orders L. For selected coupling values in the 
non-Borel-summable region, g = — .01, —.1, — 1, —10, we want to see the error as a function of the order. We want 
to find from this model system the rule for selecting systematically the best zero of P( L V) solving Eq. tfFZty . which 
leads to the optimal value of the variational parameter f2. For this purpose we plot the variational results of all zeros. 
This is shown in Fig. |2 where the logarithm of the deviations from the exact value is plotted against the order L. 
The outcome of different zeros cluster strongly near the best value. Therefore, choosing any zero out of the middle 
of the cluster is reasonable, in particular, because it does not depend on the knowledge of the exact solution, so that 
this rule may be taken over to realistic cases. 

We wish to emphasize, that for the Borel-summable domain with g > 0, variational perturbation theory has the 
usual fast convergence in this model. In fact, for g = 10, probing deeply into the strong-coupling domain, we find 
rapid convergence like A(L) ~ 0.02 exp (— 0.73L) for L — > oo, where A(L) = log \zi&) — Z oxact | is the logarithmic error 
as a function of the order L. This is shown in Fig. EI Furthermore, the strong-coupling coefficients 6; of Eq. (|2.3|l are 
reproduced quite satisfactorily. Having solved P^ L \a) = for a, we obtain Q( L \g) by solving Eq. (|2.21|) . Inserting 

this and l|2.22|l into (|2.20|) . we bring g 1 / 4 zi^(g) into a form suitable for expansion in powers of g~ 1 ^ 2 ■ The expansion 
coefficients are the strong-coupling coefficients b\ L ^ to order L. In Fig. 2] we have plotted the logarithms of their 



5 





FIG. 2: Logarithm of deviation of the variational results from exact values log \Z^ — Z cxact | plotted against the order L for 
different g < in the non-Borel-summable region. All complex optimal 57's have been used. 



absolute and relative errors over the order L, and find very good convergence, showing that variational perturbation 
theory works well for our test-model Z(g). 

A better selection of the optimal Q values comes from the following observation. The imaginary parts of the 
approximations near the singularity at g — show tiny oscillations. The exact imaginary part is known to decrease 
extremely fast, like exp (l/4<?), for g — > 0— , practically without oscillations. We can make the tiny oscillations more 
visible by taking this exponential factor out of the imaginary part. This is done in Fig. [SJ The oscillations differ 
strongly for different choices of Q( L ^ from the central region of the cluster. To each order L we see that one of 
them is smoothest in the sense that the approximation approaches the singularity most closely before oscillations 
begin. If this is chosen as the optimal one, we obtain excellent results for the entire non-Borel-summable region 
g < 0. As an example, we pick the best zero for the L = 16th order. Fig. [5] shows the normalized imaginary part 
calculated to this order, but based on different zeros from the central cluster. Curve C appears optimal. Therefore 
we select the underlying zero as our best choice at order L = 16 and calculate with it real and imaginary part for the 
non-Borel-summable region —2 < g < —.008, to be compared with the exact values. Both are shown in Fig. where 
we have again renormalized the imaginary part by the exponential factor exp (— l/4g). The agreement with the exact 
result (solid curve) is excellent as was to be expected because of the fast convergence observed in Fig. [3 It is indeed 
much better than the strong-coupling expansion to the same order, shown as a dashed curve. This is the essential 
improvement of our present theory as compared to previously known methods probing into the tunneling regime 19] . 

This non-Borel-summable regime will now be investigated for the quantum-mechanical anharmonic oscillator. 
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FIG. 3: Logarithm of deviation of variational results from exactly known value A(L) = log \Z V J — Z exact \, plotted against the 
order L for g — 10 in Borel-summable region. The real positive optimal Q have been used. There is only one real zero of the 
first derivative in every odd order L and none for even orders. There is excellent convergence A(L) ~ 0.02 exp (— 0.73L) for 
L — > co. 
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FIG. 5: Normalized imaginary part Im[.Zvar' (g) exp (— l/4g)] as a function of g based on six different complex zeros (thin 
curves). The fat curve represents the exact value, which is Z oxact (g) ~ —0.7071 + .524g — 1.78(? 2 . Oscillations of varying 
strength can be observed near g — 0. Curves A and C carry most smoothly near up to the origin. Evaluation based on either 
of them yields equally good results. We have selected the zero belonging to curve C as our best choice to this order L = 16. 
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FIG. 6: Normalized imaginary part lm[z!^P (g) exp (— 1/4<?)] to the left and the real part Re[.Zvar' (<?)] to the right, based on 
the best zero C from Fig- El are plotted against log \g\ as dots. The solid curve represents the exact function. The dashed curve 
is the 16th order of the strong-coupling expansion Z^liugid) °f equation l|2.3|l . 



III. TUNNELING REGIME OF QUANTUM-MECHANICAL ANHARMONIC OSCILLATOR 

The divergent weak-coupling perturbation expansion for the ground state energy of the anharmonic oscillator in 
the potential V(x) — x 2 /2 + g x 4 to order L 

L 

<Lk(5)=E a l9 \ (3.1) {WEAK} 

;=o 

where ai = (1/2, 3/4, —21/8, 333/16, —30885/128, ...), is non-Borel-summable for g < 0. It may be treated in 
the same way as Z{g) of the previous model, making use as before of Eqs. (|2.2U|) - (|2.23|) . provided we set p = 1 and 
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oj = 2/3, so that q = 3, accounting for the correct power behavior Eo(g) oc g 1 / 3 for g — > oo. According to the principle 
of minimal dependence and oscillations, we pick a best zero for the order L = 64 from the cluster of zeros of Pt(cr), 
and use it to calculate the logarithm of the normalized imaginary part: 

f(g) := log [^gj2 - l/3 fl . (3.2) {IM64} 

This quantity is plotted in Fig. [7| against log(— g) close to the tip of the left-hand cut for — .2 < g < —.006. Comparing 
our result to older values from semi-classical calculations |20| 

f(g)^b 1 g-b 2 g 2 + b 3 g 3 -b 4 g 4 + ... 1 (3.3) {zj} 

with 

&i = 3.95833 b 2 = 19.344 b 3 = 174.21 b 4 = 2177 , (3.4) {ZJl} 

shown in Fig. as a thin curve, we find very good agreement. This expansion contains the information on the 
fluctuations around the critical bubble. It is divergent and non-Borel-summable for g < 0. In Appendix A we have 
rederived it in a novel way which allowed us to extend and improve it considerably. 

Remarkably, our theory allows us to retrieve the first three terms of this expansion from the perturbation expansion. 
Since our result provides us with a regular approximation to the essential singularity, the fitting procedure depends 
somewhat on the interval over which we fit our curve by a power series. A compromise between a sufficiently long 
interval and the runaway of the divergent critical-bubble expansion is obtained for a lower limit g > —.0229 ± .0003 
and an upper limit g = —0.006. Fitting a polynomial to the data, we extract the following first three coefficients: 

6! = 3.9586 ± .0003 63 = 19.4 ±.12 6 3 = 135±18. (3.5) {INST} 

The agreement of these numbers with those in l|3.3|l demonstrates that our method is capable of probing deeply into 
the critical-bubble region of the coupling constant. 



Further evidence for the quality of our theory comes from a comparison with the analytically continued strong- 
coupling result plotted to order L = 22 as a fat curve in Fig. [7| This expansion was derived by a procedure of 
summing non-Borel-summable series developed in Chapter 17 of the textbook 4J. It was based on a two-step process: 
the derivation of a strong-coupling expansion of the type l|2.3[l from the divergent weak-coupling expansion, and an 
analytic continuation of the strong-coupling expansion to negative g. This method was applicable only for large 
enough coupling strength where the strong-coupling expansion converges, the so-called sliding regime. It could not 
invade into the tunneling regime at small g governed by critical bubbles, which was treated in [4| by a separate 
variational procedure. The present work fills the missing gap by extending variational perturbation theory to all g 
arbitrarily close to zero, without the need for a separate treatment of the tunneling regime. 

It is interesting to see, how the correct limit is approached as the order L increases. This is shown in Fig. [HI based 
on the optimal zero in each order. For large negative <?, even the small orders give excellent results. Close to the 
singularity the scaling factor exp (— l/3g) will always win over the perturbation results. It is surprising, however, how 
fantastically close to the singularity we can go. 



IV. DYNAMIC APPROACH TO THE CRITICAL-BUBBLE REGIME 



Regarding the computational challenges connected with the critical-bubble regime of small g < 0, it is worth to 
develop an independent method to calculate imaginary parts in the tunneling regime. For a quantum-mechanical 
system with an interaction potential g V(x), such as a the harmonic oscillator, we may study the effect of an infinites- 
imal increase in g upon the system. It induces an infinitesimal unitary transformation of the Hilbert space. The new 
Hilbert space can be made the starting point for the next infinitesimal increase in g. In this way we derive an infinite 
set of first order ordinary differential equations for the change of the energy levels and matrix elements (for details 
see Appendix B): 

K($)=Vnn{g), (4.1) {DYN} 

n y mk{g)V kn {g) sr-y V mk (g)V kn (g) 

This system of equations holds for any one-dimensional Schroedinger problem. Individual differences come from the 
initial conditions, which are the energy levels E n (0) of the unperturbed system and the matrix elements V nm (0) of the 
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FIG. 7: Logarithm of the imaginary part of the ground state energy of the anharmonic oscillator with the essential singularity 
factored out for better visualization, 1(g) — log ^y/—irg/2 E^l r (g)^ —l/3g, plotted against small negative values of the coupling 
constant —0.2 < g < —.006 where the series is non-Borel-summable. The thin curve represents the divergent expansion around 
a critical bubble of Ref. i20ll . The fat curve is the 22nd order approximation of the strong-coupling expansion, analytically 
continued to negative g in the sliding regime calculated in Chapter 17 of the textbook Q. 
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FIG. 8: Logarithm of the normalized imaginary part of the ground state energy log (^—ng/2 E a y ar (g)) — 1/3(7, plotted against 
log (— (?) for orders L = 4, 8, 16, 32 (curves). It is compared with the corresponding results for L — 64 (points). This is 
shown for small negative values of the coupling constant —0.2 < g < —.006, i.e. in the non-Borel-summable critical-bubble 
region. Fast convergence is easily recognized. Lower orders oscillate more heavily. Increasing orders allow closer approach to 
the singularity at g = 0— . 
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FIG. 9: Imaginary part of the ground state energy of the anharmonic oscillator as solution of the coupled set of differential 
equations 14. II . truncated at the energy level of n = 64 (points), compared with the corresponding quantity from the L = 64th 
order of non-Borel-summable variational perturbation theory (curve), both shown as functions of the coupling constant g. 



interaction V(x) in the unperturbed basis. For a numerical integration of the system a truncation is necessary. The 
obvious way is to restrict the Hilbert space to the manifold spanned by the lowest N eigenvectors of the unperturbed 
system. For cases like the anharmonic oscillator, which are even, with even perturbation and with only an even state 
to be investigated, we may span the Hilbert space by even basis vectors only. Our initial conditions are thus for 
n = 0, 1, 2, N/2: 



E 2n (0) =2n + l/2 
Van,2m =0 if TTi < or m > N/2 
V 2n>2 »(0) =3(8n 2 + 4n+l)/4 

^2„,2„± 2 (0) =(477 + 3)v / (277+l)(277 + 2)/2 



^2n,2n± 4 (0) = y/(2n + l)(2n + 2){2n + 3)(2n + 4)/4 



(4.3) 
(4.4) 
(4.5) 
(4.6) 

(4.7) 
(4.8) 



For the anharmonic oscillator with a V(x) = x A potential, all sums in equation (|4.1|) are finite with at most four terms 
due to the near-diagonal structure of the perturbation. 

In order to find a solution for some g < 0, we first integrate the system from to \g\, then around a semi-circle 
g = \g\ exp (iip) from tp = to ip = tt. The imaginary part of E$(g) obtained from a calculation with TV = 64 is shown 
in Fig. refVIx, where it is compared with the variational result for L = 64. The agreement is excellent. It must 
be noted, however, that the necessary truncation of the system of differential equations introduces an error, which 
cannot be made arbitrarily small by increasing the truncation limit N . The approximations are asymptotic sharing 
this property with the original weak-coupling series. Its divergence is, however, reduced considerably, which is the 
reason why we obtain accurate results for the critical-bubble regime, where the weak-coupling series fails completely 
to reproduce the imaginary part. 



{INI} 



APPENDIX A 



We determine the ground state energy function E$(g) for the anharmonic oscillator on the cut, i.e. for g < in the 
bubble region, from the weak coupling coefficients a/ of equation (|3.1|l . The behavior of the a/ for large I can be cast 
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into the form 

L 

ai/oi_i r J . (5.1) 

i=-l 

The /3j can be determined by a high precision fit to the data in the large I region of 250 < I < 300 to be 

( 3 95 113 391691 40783 1915121357 10158832895 70884236139235 
P~ x ' °' h ~\ ' ~2' 24' ~' 3456 ' 48 ' 248832 ' 124416 ' 71663616 ' ^ ' ' 

60128283463321 286443690892 144343264152266 351954117229 2627843837757582 



4478976 ' 1423 43743 ' 6 2339 

230619387597863 12122186977970425 41831507430222441029 



10 24 3550 

where the rational numbers up to j — 6 are found to be exact, whereas the higher ones are approximations. 
Equation l|5.1[) can be read as recurrence relation for the coefficients a;. Now we construct an ordinary differential 
equation for E(g) := Eq^ a j,(<?) from this recurrence relation and find: 



d N L L+X 



dg 



L+i f d V 



9 

3=0 



E(g) = . (5.3) 



All coefficients being real, real and imaginary part of E(g) each have to satisfy this equation separately. The point 
g = 0, however, is not a regular point. We are looking for a solution, which is finite when approaching it along the 
negative real axis. Asymptotically E{g) has to satisfy E(g) ~ exp(l/g ( 9_ 1 ) = exp(l/3g). Therefore we solve l|5.3ll 
with the ansatz 



E{g) = g a exp^-X>(- S )^ 



(5.4) 



to obtain a = —1/2 and 



. 95 619 200689 2229541 104587909 7776055955 9339313153349 172713593813181 
i,2,3,... -i 24> I2" > 1152 ' 1024 ' 3072 ' 12288 ' 688128 ' 524288 ' ^ ' ' 
1248602386820060039 14531808399402704160316631 12579836720279641736960567921 



139886592 54391637278720 1435939224158208 

109051824717547897884794645746723 45574017678173074497482074500364087 



348951880031797248 ' 3780312033677803520 

This is in agreement with equation 13.4fl and an improvement compared to the WKB results of [2(|. Again, the first 
six rational numbers are exact, followed by approximate ones. 

VI. APPENDIX B 

Given a one-dimensional quantum system 

(H +gV)\n,g)=E n (g)\n,g) (6.1) 

with Hamiltonian H = Hq + g V, eigenvalues E n (g) and eigenstates \n, g) we consider an infinitesimal increase dg in 
the coupling constant g. The eigenvectors will undergo a small change: 

\n,g + dg) = \n,g) + dg ^u nk \k,g) (6.2) 



so that 

d_ 

dg 



= u ™k\k,g) ■ (6.3) 
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Given this, we take the derivative of (|6.1(l with respect to g and multiply by (m,g\ from the left to obtain: 

(m,g\V-E' n (g)\n,g) = Y / u nk (m,g\H Q + gV-E n (g)\k,g) . (6.4) 

Setting now m — n and m ^ n in turn, we find: 

Kid) =Vnn(g) (6.5) 
V mn (g) =u nm (E m (g) - E n (g)) , (6.6) 

where V mn (g) = (m,g\V\n,g). 

Equation (|6.5|) governs the behavior of the eigenvalues as functions of the coupling constant g. In order to have a 
complete system of differential equations, we must also determine how the V mn {g) change, when g changes. With the 
help of equations (|6.3|l and (|6.6|) , we obtain: 

VLn = u *mk( k :9\V\ n ,g) + u nk{m, g\V\k, g) (6.7) 

k^m k^'n 

Equations l|6.5|l and (|6.8|) together describe a complete set of differential equations for the energy eigenvalues E n (g) 
and the matrix-elements V nm (g). The latter determine via l|6.6|l the expansion coefficients u mn {g). Initial conditions 
are given by the eigenvalues E n (0) and the matrix elements V nm (0) of the unperturbed system. 
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